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£■""». , Abstract 

o ■ 

{Nj ■ The basic concepts and equations of classical fluid mechanics are presented in the 

form necessary for the formulation of Newtonian cosmology and for derivation and 
\ analysis of a system of the averaged Navier-Stokes-Poisson equations. A special 

""^ ■ attention is paid to the analytic formulation of the definitions and equations of 

Q" 1 moving fluids and to their physical content. 

U '. 

> : 1 Introduction 

In looking for adequate gravitational field equations applicable in the cosmological settings 
! l] there has been a number of attempts to average Einstein equations (for a review and 
references see [fj, and recent papers ||, 0). Despite four decades of research in 
this field since the pioneering paper ||, there is still no definite consensus on the results 
and approaches to the averaging problem in general relativity and cosmology. It does 
reflect to a certain extent the fundamental character of this problem. The averaging 
problem in general relativity has actually originated in cosmology |T]]-[§], 0-[§|, but the 
main questions remain unanswered. A definite and unambiguous picture of the averaged 
dynamics of the universe is not available as yet. 

In view of the difficulties mentioned above, there have been some attempts to con- 
struct an averaged theory of the Newtonian self-gravitating fluid to describe on its basis 
various aspects of dynamics of homogeneous and inhomogeneous Newtonian cosmologies, 
such as backreaction, large-scale structure formation, scaling properties (for a review and 
references see M, [ITU] and recent papers |TTJ-|TJ§). Studying the Newtonian gravity and 



cosmology is of importance for three main reasons: 
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(i) The physics of Newtonian self-gravitating configurations provides an adequate the- 
oretical framework for modelling the structure formation, galaxy dynamics, clustering 



processes, etc. in cosmology |L4 



(ii) It is well-known that Newtonian cosmology has very similar physics in many as- 



pects with the general relativistic cosmology [15], [16] 



(iii) Newtonian problems are simpler to deal with than the general relativistic ones and, 



given similarities and differences between them [17 , , one can gain important intuition 
and experience to approach the similar general-relativistic cosmo logical problems.. 

In order to derive a system of the averaged Navier-Stokes-Poisson equations and to 
make use of the system for analysis of inhomogeneous Newtonian universes, it is signifi- 
cant to present the basic concepts and equations of classical fluid mechanics in the form 
necessary to approach this problem. In this presentation a special attention will be paid 
to the analytic formulation of the definitions and equations of moving fluids and their 
physical content. 

The structure of this paper is as follows. Chapter 2 gives a brief account on the 
interconnection between fluid mechanics and Newtonian cosmology. In Chapter 3 some 
approaches to averaging inhomogeneous Newtonian cosmologies are reviewed and dis- 
cussed. The main concepts regarding the fluid motion and the Eulerian and Lagrangian 
pictures of moving fluids are presented in Chapters 4. Chapter 5 presents the field de- 
scription of moving fluids which provides a framework for an analytic formulation of fluid 
motion as the theory of the integrals of partial differential equations. The Euler expan- 
sion formula for the change of the volume of an infinitesimal element of fluid during its 
motion is proved in Chapter 6. The Reynolds transport theorem following from the Euler 
expansion formula for the change of an integral of a fluid quantity over an arbitrary fluid 
volume moving along the fluid paths is proved in Chapter 7. In Chapter 8 the fluids 
possessing the property of inertia are introduced and the equation of continuity for the 
fluid density is derived as a consequence of the fluid mass conservation. The Principle of 
conservation of linear momentum is formulated in Chapter 9 and the Cauchy equation of 
motion is obtained as a consequence of this conservation. The structure of the system of 
the Cauchy equation of motion and the equation of continuity for the description of fluid 
motion is also analyzed here. The kinematic quantities to characterize the change of a 
fluid element during its motion, namely, the expansion scalar, the shear tensor and the 
vorticity vector, are introduced and discussed in Chapter 10. Chapter 11 is devoted to 
analysis of the conservation of energy for fluids and formulation of the laws of thermody- 
namics for fluids with a special attention paid to the following energy transfer equations 
and the definition of the thermodynamical variables for the compressible and incompress- 
ible fluids. The next Chapter gives a discussion of the conservation of angular momentum 
and its relation with the structure of the fluid stress tensor. Chapter 13 provides a list of 
different forms of the Cauchy equation of motion in dependence on the structure of fluids. 
In particular, the Euler equation for the perfect fluid, the Stokes equation for the Stokean 
fluid and the Navier-Stokes equation of classical hydrodynamics for the Newtonian fluid 
are presented and discussed. The physical status of the Navier-Stokes equation as the 
macroscopic equation governing the dynamics of gases and liquids is also discussed. 
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Conventions and notations are as follows. All functions / = f(x l ,t) are defined on 3- 
dimensional Euclidean space E 3 in the Cartesian coordinates {x 1 } with Latin space indices 
i,j, k, ... running from 1 to 3, and t is the time variable. The Levi-Civita symbol is 
defined as £123 = +1 and e 123 = +1 and b % \ 8j and Sij are the Kronecker symbols. The 
symmetrization of indices of a tensor T l - k = Tj k (x l , t) is denoted by round brackets, 7?.^ = 
\ \Tjk + ^fcj); an d the antisymmetrization by square brackets, T?. fc , = ~ (Tj fc — T^-). A 
partial derivative of Tj fc with respect to a spatial coordinate or time t is denoted either 
by comma, or by the standard calculus notation, T l - kt = dT l - k /dt and T l - k ^ = dT l - k /dx l , 
the fluid velocity is u % = u l (x\t), and the material (total) derivative is T l - k = dTL/dt = 
dT; k /dt + u l dT; k /dx l = 7},, + „%,,. 

2 Fluid Mechanics and Newtonian Cosmology 

The dynamics of a compressible Newtonian fluid evolving in its own Newtonian gravita- 
tional potential field is known to be governed by the Navier-Stokes-Poisson system (see 
19 1, |20| and plj| , [p^j for the system of equations in the framework of the so-called 



3+1 covariant formulation of general relativistic cosmology and references therein). It 
consists of the equation of continuity which states the conservation of mass of a fluid 
element, the Navier-Stokes equation which is the Newton second law for the element of 
fluid stating the conservation of the fluid element momentum, the Poisson equation for 
the fluid gravitational potential and an equation of state of the fluid. The fluid flow is 
described by its velocity, density, pressure, or temperature, and the molecular transport 
coefficients of dynamic and bulk viscosity and/or thermal conductivity. Generally, the 
evolution of a fluid element can be characterized by change of its kinematic characteris- 
tics: the expansion scalar of the fluid element characterizes the change of its volume, the 
shear tensor determines the distortion of the fluid element leaving its volume unchanged 
and the vorticity vector represents its rigid rotation. The Navier-Stokes equation and 
the equation of continuity can be therefore rewritten explicitly as a system of evolution 
equations for the expansion, the shear, the vorticity and the fluid density. In case of New- 
tonian cosmology the system of the Navier-Stokes-Poisson equations for a compressible 
Newtonian fluid without viscosity and heat transfer is assumed to describe the Newtonian 
universes, that is, possibly spatially infinite distributions of the evolving self-gravitating 
cosmological fluid governed by the Navier-Stokes-Poisson equations. For homogeneous 
and isotropic Newtonian universes only time-dependent functions of the expansion, the 
pressure and the density remain essential. The evolution equation for the expansion scalar 
then describes the dynamics of evolving homogeneous and isotropic Newtonian universes. 
This equation is known to have the same form as the Raychaudhuri equation in general 



relativity |[23||, [24] describing the evolution of the expansion scalar of a 4-dimensional flow 



of the cosmological fluid in homogeneous and isotropic general relativistic cosmology. It 



is equivalent to the Friedmann equation [25] for a scale factor of an expanding universe in 
both Newtonian and general relativistic cases. Therefore, studying Newtonian cosmology 
proved its importance for understanding some aspects of general relativistic cosmology 
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and revealed remarkable similarities between both theoretical approaches to studying the 
cosmological models of the real Universe 



15 , 16 



In case of inhomogeneous and anisotropic Newtonian universes, when the whole system 
of the equations must be considered, one can study the dynamics of a Newtonian universe 
which is smooth, homogeneous and/or isotropic, on large scales. In the general relativis- 
tic setting an inhomogeneous cosmological model which has the Friedmann-Lemaitre- 
Robertson- Walker (FLRW) space-time geometry [p5fl-[29[j on average is generally believed 
to be more realistic than the standard FLRW model of our Universe M. An essential 
dynamical factor which must affect the character of dynamics of such an inhomogeneous 
universe is that it takes into account the effect of the backreaction of inhomogeneities on 
the global expansion of the universe (for a review and references see and recent papers 
12 1, [fl3 1 , [[31~| ). In case of Newtonian cosmology a similar setting is also of great interest 



for understanding the cosmological backreaction due to correlations of a self-gravitating 
cosmological fluid. To construct such a model a suitable averaging procedure must be 
applied to the Navier-Stokes-Poisson system, to bring about a set of equations for the 
averaged expansion, shear, vorticity, pressure, fluid density and gravitational potential 
together with their correlation functions involved. 

The nonrelativistic kinematics and dynamics of fluid flows have been intensively stud- 
ied since the eighteenth century (for a review and references see [[32| , [[33|] ) , which resulted 
in discovering the Navier-Stokes equation as the fundamental equation of classical hydro- 
dynamics H34 
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This equation successfully describes the motion of various fluids and 
gases and other transport phenomena under a wide class of physical conditions including 
viscosity and heat transfer ||20|| . The existence of two different types of flow, called lam- 
inar and turbulent, has been established in the middle of the nineteenth century. The 
pioneering systematic investigations of turbulence by Rayleigh and Reynolds in the last 
decades of the nineteenth century led to understanding that the correlation functions of 
the fluid velocity and the fluid density must be recognized as the fundamental charac- 
teristics of turbulence (for a review and references see [ID], [3~6"|-[3~S|). The first system 
of differential equations for the correlation functions containing the single-point first and 
second moments of the fluid fields has been derived by Reynolds |$9| using the time and 
space averagings of the system of the continuity and the Navier-Stokes equations for the 
case of incompressible fluids. An analytical formulation of the problem of turbulence 
has been given first by Keller and Friedmann |^0[ who developed a general method for 
obtaining the differential equations for the moments of arbitrary order by using the sys- 
tem of the continuity and the Navier-Stokes equations. The Reynolds equations are the 
first equations in the total infinite system of the Keller-Friedmann equations. Analysis 
of these systems of equations has resulted in remarkable progress in our understanding 
of the phenomenon of turbulence to establish the statistical nature of turbulence, to de- 
velop the theory of homogeneous and isotropic turbulence, to propose the hypothesis of 
locally isotropic turbulence and developed turbulence 



0, 0-|8|. However, despite 
the significant developments the problem of turbulence remains essentially unsolved since 
even such fundamental questions as which flow should be considered turbulent and which 
should not, what is a mechanism of the turbulent mixing, how the developed turbulence 
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occurs, etc. do not have satisfactory explanation as yet. The difficulties originate partly 
in the mathematics of the Navier-Stokes and Reynolds equations where the questions of 
existence, uniqueness, and smoothness of solutions are not mainly resolved On the 
other hand, theoretical approaches to formulate adequate physical models of real physical 
phenomena exhibiting turbulent processes are often very limited due to the absence of 
sufficient experimental evidence about a turbulent regime of interest W2 . 



3 Approaches to Averaging out Inhomogeneous Uni- 
verses 

As it has been pointed out above, the physics of Newtonian self-gravitating fluid con- 
figurations is widely employed for modelling cosmological gravitating structures and the 
Universe a whole. The idea of Chaos, or Cosmic Turbulence, as the origin of our Universe 
goes back to Kepler, Descartes and Laplace . Development of the hot Big Bang Model 



and attempts to understand the origin of galaxies has led to invoking the hydrodynamic 
turbulence theory in modelling the galaxy formation ||44|| , p5[ as opposed to the idea of 
the origin of galaxies due to the gravitational instability (see [fL4H , ||46f for discussion and 



references). It has been shown [47|, in particular, that the physical conditions in the early 
universe could allow the Primeval Turbulence regime and the appearance of eddies to 
produce galaxies later on. 

However, application of the Navier- Stokes- Poisson system and its suitably averaged 
version for analysis of the evolution of inhomogeneous and anisotropic Newtonian uni- 
verses has not been accomplished in any considerable depth. Such important questions 
as backreaction, large-scale structure formation, scaling properties, dynamics of universe 
which is homogeneous and isotropic on average and other dynamical issues related to the 
character of the cosmological evolution when the gravitational correlation effects are taken 
into account are still open. In this setting the system of equations for the 3+1 covariant 
formulation of general relativistic cosmology has been hardly studied either. 

A study of the averaged-over-ensemble self-gravitating Newtonian perfect fluid in the 
Eulerian formulation has been carried out first in |§ under assumption of the homogeneous 



isotropic turbulence Despite this strong requirement which leads to a considerable 
simplification of the system of the averaged Navier- Stokes- Poisson equations, the New- 
tonian cosmological models exhibit nontrivial dynamics. The evolution equation for the 
mean-square correlation function of vorticity has been derived and its solution has been 
shown to have a nontrivial behavior to blow up in a finite time for an expanding back- 
ground in case of domination of the fragmentation of eddies (inhomogeneities), and to 
blow up always in a finite time for a contracting background. Thus, in dependence of a 
value of the initial vorticity the averaged evolution of an initially inhomogeneous Newto- 
nian cosmological model can develop either a turbulence mixing regime, or a regime of 
the natural decay of inhomogeneities. 

Analysis of the general relativistic Raychaudhuri equation has been attempted in JK| 
for the case of closed universes which are globally splitted into a topological product of 
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compact, orientable, differentiable 3-manifolds and the real line which represents time. 
These 3-manifolds can be naturally represented as a family of disjoint spacelike hyper- 
surfaces parametrized by time, and there exists a future-oriented normal unit 4-vector 
field tangent to a normal congruence of curves which are orthogonal trajectories of the 
hypersurface family. A general relativistic analogue of the Reynolds transport theorem 



of fluid mechanics [|32|| , [53], \W\ for the rate of change of the integral of the expansion 
scalar of this 4-vector flow over a hypersurface moving with the flow has been derived. 
A consequent formula for the rate of change of the ratio of this integral to the volume of 
the corresponding hypersurface, which has been called the "average expansion" , was used 
to find out if this "average expansion" has a zero volume singularity characteristic to the 
FLRW models. As a result of this analysis it has been suggested that such a singularly 
can be avoided at all if during some epochs the universe is expanding highly nonuniformly. 

This result, however, should not be misinterpreted. Since by its kinematic meaning a 
transport formula for a quantity is essentially an integral formulation of the corresponding 
differential law for this quantity, there is no new either kinematic, or dynamical informa- 
tion which can be gained by using such an integral formulation. So, for example, the 
conservation of the mass of fluid in a material volume moving with the fluid leads to 
the equation of continuity which is well-known to be a necessary and sufficient condition 
for a motion to conserve the mass of an arbitrary moving volume ||32|| . On the contrary, 
the averaged continuity equation obtained by the application of a properly defined space, 
time or ensemble average and the corresponding formulae for interchanging of operations 
of taking a derivative and an averaging is known to be a new differential dynamical equa- 



tion, the continuity equation for the turbulent flow JT9|, p8[ . Therefore, the conclusions 
derived in ||10| on the basis of using a ratio of the integral of the expansion scalar over 
an arbitrary hypersurface moving along the normal 4-vector flow to the volume of this 
region as a space average of the expansion and using the transport formula for this ratio 
as a formula for interchanging of the time derivative and this "averaging" cannot be con- 
sidered as a consistent averaging procedure. There is no new information in the derived 
integral formula as compared with the corresponding local differential relation which is 
the general relativistic Raychaudhuri equation. 

Some recent approaches to construct an averaged theory of self-gravitating Newtonian 
fluid to study the dynamical properties of the homogeneous and inhomogeneous Newto- 
nian universes, such as backreaction, large-scale structure formation, scaling properties 
(see |ll|]-[[T3| and references therein), have also made use of the similar ratio described 



above pi as a "space average" of the expansion scalar and other kinematic quantities 



for the case of 3-dimensional cosmological fluid. A nonrelativistic version of the transport 



formula as in [|Hj] which is known from the Reynolds transport theorem ||32|| , [p3fl , |49 
modified by division by the region volume has been taken as a formula for interchanging 
of the time derivative and this "averaging" and a set of the transport relations for the 
kinematic quantities has been derived. This set of transport relations was regarded as the 
averaged equations for the kinematic quantities. In particular, the transport equation for 
the expansion scalar has been taken as the averaged Newtonian version of Raychaudhuri 
equation to analyze the character of the averaged dynamics. Again, such a procedure of 
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space averaging is not proper and all conclusions are achievable by analyzing the initial 
system of the differential equations for kinematic quantities. 



4 The Eulerian and Lagrangian Pictures of Fluid Mo- 
tion 

The material fluid is considered to be a continuous physical substance distributed over 
a region of the 3-dimensional Euclidean space, which is assumed to possess a positive 
volume. A fluid flow, or a fluid motion, is represented mathematically by a one-parameter 
family of mappings of a 3-dimensional region filled with the fluid into itself, the mapping 
being usually assumed to be 1-1 and at least thrice continuously different iable in all 
variables |32fl ) [ j33|| , |50| , |5l| except possibly at certain singular surfaces, curves or points 



where a separate analysis is required. The real parameter t describing the mappings of 
fluid regions is identified with the time with its domain supposed to be — oo < t < +00 
where t = is an arbitrary initial instant. 

In order to describe the fluid motion analytically let us introduce a fixed rectangular 
coordinate system {x 1 }, each coordinate triple of which is referred to as a position of 
an infinitesimal element of a fluid distribution, or configuration, often called a "fluid 
particle" . If one observes such a fluid particle moving with the fluid, it is easy to detect 
that the particle which was in the position has moved to position {x 1 } at a later 
time. Without loss of generality, one can take the first instance to be t = and if the 
later time is t the fluid flow is represented by the transformation 

x* = x\e,t). a) 

If £ l is fixed while t varies, Eq. ([I]) specifies the path of the fluid particle 

X * = x\t) = x*(e,t)ie =const , x*(o) = e, (2) 

initially placed at £\ If time t is fixed, Eq. ([!]) determines a transformation 

X* = X\e) = ^(e',%=const, Z*(0) = C, (3) 

of the region initially occupied by the fluid into its position at time t. The transformation 
(|IJ) is precisely the above-mentioned one-parameter family of mappings defining the fluid 
motion 

x i :x i (0) = ?^x i (?,t). (4) 

Since functions x l (^,t) are single- valued and thrice differentiable in all variables, the 
Jacobian of the transformation (JI|) does not vanish, 

J= S^ =det fggVo, (5) 
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and there exists a set of inverse functions 

e = e(x\t) (6) 

to define the initial positions of the fluid particle which is at any position {x 1 } at a moment 
of time t. This inverse transformation is also 1-1 and thrice different iable in all variables. 
Due to Eq. @ a fluid flow can be equivalently described by a set of initial positions of 
its particles in dependence of their positions in later times. 

Thus, there are two equivalent ways to describe the fluid motion. 

The Eulerian description by the spatial variables (x\t) This picture presents the 
state of motion at each position {x 1 } such that throughout all time a given set of coordi- 
nates {x 1 } remains attached to a fixed position, serving as an identification for it - the 
spatial variables are called the Eulerian or spatial coordinates. 

The Lagrangian description by the material variables (£\t) This picture of the 
motion chronicles the history of each fluid particle such that throughout all time 
a given set of coordinates remains attached to a single fluid particle, serving to 
uniquely identify it - the material variable are called the Lagrangian, material or convected 
coordinates []• 

It should be pointed out here that the Eulerian and Lagrangian descriptions have quite 
different physical interpretations when one would like to describe a moving fluid from the 
point of view of an observer carrying out measurements by a set of physical instruments. 
In the Eulerian picture an observer is always located in a particular position {x 1 } at a 
moment of time t and such an observer views moving fluid particles passing through and 
by the position. In the Lagrangian picture an observer is always moving with a fluid 
particle which was initially at a position and such an observer views the relative 
changes in the fluid comoving with the observer's particle. 

5 The Field Description of Fluid Motion 

Although the fluid flow is completely determined either by the transformation ([]]) in the 
Eulerian description, or its inverse transformation (|6|) in the Lagrangian description, it is 
important also to be able to describe the state of motion at a given position during the 
course of time. To achieve this goal one should characterize the motion of fluid in terms of 
objects defined at a position {x 1 } at a time t. A simplest object to be considered for this 
role is the instant fluid velocity u l , functions of which can be defined in both descriptions, 
u\x\t) and u\&,t). 

Let us consider a quantity / which characterizes a fluid in the framework of the 
Eulerian or Lagrangian descriptions. Generally, any such a quantity / is a function of 
spatial variables / = f(x\t) and it is also a function of material variables / = f{£ l ,t) 

1 Though the descriptions of the fluid motion by the spatial and material coordinates are usually 
called Eulerian and Lagrangian, correspondingly, both of them are known to be due to Euler p2| , [p0| . 
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due to 

f = f[x l (e,t),t] or f = f[C{x j ,t),t] } (7) 

the functions (|7|) being related by Eqs. (|J) and (|j) considered as the laws of the change 
of the variables. Geometrically, / = f(x\t) is the value of / felt by the fluid particles 
instantaneously located at the position {x 1 } and / = t) is the value of / experienced 
at time t by the fluid particle initially located at the position Then the change in 

quantity / when the fluid moves can be characterized by the two different time derivatives 

% = and f^mA . (8) 

Ot Ot |»<=conat at Ot |£i=const 

The partial derivative df/dt gives the rate of change of / with respect to a fixed position 
{x 1 }, while the material, or convective, derivative df/dt measures the rate of change of / 
with respect to a moving fluid particle. 

The material derivative of the position of a fluid particle defined by the second Eq. 
(H) for / = x l is called its velocity u l 

AC.t) = d 4 = dj ^ ■ P) 

Ul at |g l =const 

The notion of the fluid velocity plays a fundamental role in the field formulation of fluid 
mechanics. Indeed, in most physical settings, when the actual motion of a fluid particle, 
that is either x % = x l (^,t) or £ J = ^(x 1 , t), is not given explicitly, it is more convenient 
to deal with the Eulerian form of the fluid velocity, u l = u l (x\t). Given a vector field 
of fluid velocity u % = u*(rr J ",t), which can be actually measured in laboratory, defined in 
a fluid 3-region U for a time interval At, one can determine the functions x % = x l (^,t) 
by solving the initial value problem for the system of the first order ordinary differential 
equations 

^L = «V.*)> * , (o) = f. (io) 

A unique solution x % = x t (^,t) to Eqs. ( |10D is known (see, for instance, |f52| , f53fl ) to 



exist always on U x At for at least once continuously different iable in all variables and 
bounded functions u % = u l [x 3 ,t), and it depends continuously on the initial values 
to completely characterize the fluid motion (|1|). Therefore, the description of fluid motion 
in terms of the fluid velocity is equivalent to the Eulerian or Lagrangian pictures: 

The field description by the fluid velocity in the spatial variables u l {x\t) This 
picture presents the state of a fluid in motion by its velocity field specifying the velocity 
u l = u l (x\t) of a fluid particle at each position {x 1 } at each time instance t, such that 
throughout all time a given set of velocity components {«*} remains attached to a fixed 
position, serving as an identification for a single fluid particle in dependence of its initial 
position - the velocity variables are called the field variables^ 

2 Euler and d'Alembert are considered to be the first to discover the field description of moving fluids. 
However, it was Euler who put forward the idea of formulating a field theory of fluid motion as the theory 
of the integrals of partial differential equations |B2j , [M . 
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The physical interpretation of the field description of fluid motion represents an ob- 
server as located in a particular position {x 1 } at a moment of time t to measure the 
velocities of the fluid particles moving by the observer and which were initially at {£*}■ 

The field description of moving fluid has a fundamental significance for the proper 
analytic treatment of processes occurring in physical continuous media under a wide 
class of physical conditions. It is the field description that provides the possibility of an 
analytic formulation of fluid kinematics and dynamics by finding a set of partial differential 
equations governing the fluid motion. 

With the definitions of the time derivatives (§) one can calculate the material derivative 
of any quantity f(x^,t) as 



dj_ 
dt 



0_ 

ot 

d_ 
dt 



f(C,t) 

|^ 1 =const 

f[x\e,t),t] 



§ l =const 



df 

dt\ 



+ 



=const 



df_ 

dx' 



dt 



(11] 



|^ I =const 



u l {x^t) defined 



that can be written in the following form by using the fluid velocity u l 
by Eq. (0): 

df df { df 

Jt = Tt +U d^- (12) 
This formula relates the material and spatial derivatives through the velocity field in the 
field picture of fluid motion and it expresses the rate of change in f{x\ t) with respect to 
a moving fluid particle located at a position {x 1 } at time t. An important application of 
the material derivative (O) is the definition of the fluid acceleration a 1 = a*(x J , t), 



du 



du l j du l 
dt dt dxi ' 



(13) 



as the rate of change of the velocity of fluid particle experienced during its motion. Due 
to (p~3|) the fluid acceleration a 1 is represented only in terms of the fluid velocity u l . 



6 The Euler Expansion Formula 

The Jacobian J (|5|) of the transformation (|l|) represents the dilatation of the infinitesimal 
volume of a fluid particle as it moves along its path @). Indeed, considering ([l]) as the 
coordinate transformation from the initial coordinates at time t = of a fluid particle 
to the particle coordinates {x 1 } at later time t, the change of the particle's infinitesimal 
volume is given by dV = JdVo where dV = dx 1 dx 2 dx 3 is the particle volume at the 
position {x 1 } and dVo = d^d^dl; 3 is the initial particle volume. Because the fluid motion 
is continuous and ([!]) is continuously differentiable and therefore invertible, the fluid 
dilatation, or expansion, dV/dVo, is positive and bounded 

dV , x 

J =—, 0<J<oo. (14) 
avo 
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Since the infinitesimal volume of a fluid particle undergoes, in general, changes and dis- 
tortions during its motion, the rate of change of the dilatation characterizes such changes. 
The formula for the dilatation change had been found by Euler [54| (see p2| , J33] , jSOl , 
|5~T1 for discussion). 

Theorem 1 (The Euler expansion formula) The rate of change of the fluid dilata- 
tion J ( 1?^ ) following a fluid particle path is given by 

dJ T du l d du l . . 

— = J-— , or — mj = — - . 15) 

dt dx l dt dx l v ' 

Proof. A derivation of the Euler formula is straightforward by differentiation of the 
Jacobian determinant, dJ = JdC, l /dx J d (dx l /d£-> ) , 

dJ d^ d ( dx % \ dtf du % dt? du l dx k jdu 1 . . 

QED 

Since the Euler expansion formula ( |i5|) establishes the divergence of the fluid velocity 
as a measure of the fluid dilatation, or expansion, the former is also called the fluid ex- 
pansion and denoted as 9 = du % jdx % . The Euler formula QT5D emphasizes the fundamental 
role of the fluid velocity in the field description of fluid motion. When a fluid is moving 
without change in volume for all fluid particles at all position {x 1 } throughout all time t, 

dJ du l , 

m =0 * a? = ' < 17) 

such a fluid is called incompressible. If a fluid is not incompressible, 

dJ , du % , , , 

is called compressible. 



7 The Transport Theorem of Reynolds 

An important kinematical theorem follows from the Euler expansion formula (|I5D. It is 
called the Reynolds transport theorem ||5 and it concerns the rate of change not an 
infinitesimal volume of a fluid particle, but any integral of any quantity / which is a 
function of spatial variables / = f(x\t) over an arbitrary closed fluid region, or volume, 
£(£) moving with the fluid H32| , |33[| , |]50fl , Let us consider a volume integral of f(x\ t) 



over E(t), 

F(t)= I f(x\t)dV, (19) 



E(t) 

which is a well-defined function F = F(t) of time. 
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Theorem 2 (The Reynolds transport theorem) The rate of change of the volume 
integral ( |7^ following the fluid particle paths is given by 

E(t) E(t) 

Proof. It is straightforward by transformation of ([H]) into the integral over the fixed 
initial volume So = S(0) in the material coordinates 

Jt F{t) = Jt J /W,*),*]^o= / {^/W.^.tR + ZW.^tl^Wo (21) 

So So 

which immediately brings ( p0| ) by applying the Euler expansion formula (pT5|). QED 

By substitution (|12|) for the material derivative df/dt in the integrand of (p0[), one can 
obtain the Reynolds transport theorem in the following form 



dl 



(22) 



£(t) E(t) 



There is yet another useful form of the last formula obtained by application of the Gauss 
theorem 

i S fdV = I d i iv + I lutiSk (23) 

£(t) E(t) 8E(t) 

where the 2-surface element is a vector equal to the infinitesimal 2-square and normal 
to the surface a(t) = <9£(t) enclosing the volume The formula (^3|) admits an 

immediate physical interpretation as expressing the rate of change of a quantity / within 
a moving fluid volume as the integral of the rate of change in this quantity itself at a 
point plus the net flow of / over the surface bounding the volume. 

One particular case of the Reynolds transport theorem is of importance. For the 
constant quantity / = const the Reynolds transport theorem (|2"0|) gives a formula for the 
change of value of the volume Vs(t) of an arbitrary moving fluid region £(£), 



Vfc(t) = J dV. (24) 
s(t) 

Corollary 1 (The change of volume) The rate of change of the value of the volume 
Vs(i) P^ j of an arbitrary moving fluid region £(i) is determined by the formula 

d V^(t)= f^dV= [ u k dS k . (25) 



dt J dx k 

E(t) 9E(t) 
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On the basis of the Euler formulas [ 50 1 , |55| for the transformation of the infinitesimal 
line element of a fluid particle path between a spatial line element dx % and a material line 
element d£ l , 

<% k lLA?,t) and d C = dx k ^C(x^t), (26) 



dx 1 



and the Nanson formulas 



01, 156 



for the transformation of the infinitesimal surface 
elements for a 2-dimensional surface made of fluid particle paths between a spatial surface 
element dSi and a material surface element dS^, 



ds k ^x k (e,t) 



d 

JdS (0)i and dS i0)k —^ k (x 3 ,t) = JdSi, 



(27) 



one can derive the corresponding line and surface integral transport relations [33|, [50 



Slj similar to the Reynolds transport theorem. 

The Reynolds transport theorem is an integral statement essentially equivalent to the 

Euler expansion formula. The same true for the corresponding line and surface integral 

theorems and the Euler and the Nanson formulae (^), (PD- 
As it has been discussed in Introduction, the Reynolds transport theorem has been 

recently utilized for a definition of a space averaging over fluid volumes. This issue will 

be discussed in |57 



8 The Equation of Continuity 

The definition of the fluid velocity (|9|), the system of partial differential equations (|1(]) for 
determination of the equation of motion of fluid x % = x l {^ , t), the Euler expansion formula 
( |T5D and the Reynolds transport theorem (|20J ) constitute the kinematics of moving fluid. 
The formulation of kinematics does not require any hypotheses about the structure and the 
nature of fluids as physical substances. To establish the differential equations governing 
the fluid dynamics, one must formulate such physical axioms on the basis of experimental 
evidence. 

As it has been mentioned above, in the framework of the field description one must 
characterize the motion of fluids in terms of objects defined at a position {x*} at a time 
t. The fluid velocity u L has been shown to be the simplest quantity to formulate the 
kinematics of fluid. Let us now further consider fluids which possess the property of 
inertia. Such a fluid can be endowed with density function p = p(x l ,t) of the physical 
dimension [mass/volume], which enables one to define the mass M(t) of the fluid held 
within the region S(t) 

M(t) = J pdV. (28) 

E(t) 

The fluid mass ( ^8|) for any region S(t) is assumed to be always positive that means 
positivity of the density function, 

p{x\t)>0, (29) 
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everywhere. The following physical principle is taken to hold for the fluid. 

The Principle of conservation of mass The mass of fluid M(t) (|28| ) in an arbitrary 
region £(£) does not change following the fluid particle paths 

d M{t) = or ^ [ pdV = 0. (30) 



dt dt 

E(t) 



The Principle of conservation of mass immediately brings about a differential relation 
for the density function p(x%t) that is called the equation of continuity |32|], [Q, ]50 
]5T | discovered first by Euler 



Theorem 3 (The equation of continuity) The equation of continuity 

dp dv? , , 

i + ^ = ° (31 » 

is necessary and sufficient for the mass of any region following the fluid particle paths to 
be conserved (Wl 



Proof. Using the Reynolds transport theorem ( P0|) for the fluid mass M(t) (g^) gives 
the rate of change of the mass as 



dp du 

-7L + P 



dt dx k 



dV. (32) 



E(t) 



Necessity follows immediately upon substitution of the equation of continuity (|3T| ) into 
(j32|) . Sufficiency follows from the mass conservation (|30"D and (|32D by taking into account 
that the fluid region E(t) is arbitrary. QED 

Using the formula for the material derivative ([12]) gives another form of the equation 
of continuity 

| + !j(^) = °. (33) 

An important formula for the rate of change of an integral of a product of a quantity 
f(x\t) by the density function p(x\t) 



[pf](t)= j p(x\t)f(x\t)dV. (34) 

E(t) 

is due to the Reynolds transport theorem (|2Tj| ) and the equation of continuity (|?1"|). 

Corollary 2 For any quantity f(x l ,t) the rate of change of the volume integral of the 
product p(x l , t) f (x^ , t) ( |5^[ ) following fluid particle paths is given by 

JtS pSdv= S pd i dv - (35 » 

s(t) E(t) 
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The equation of continuity holds for any fluid satisfying the Principle of conservation 
of mass (p0|) . It should be pointed out here that though the fluid mass is not a kinematic 
quantity the derivation of the equation of continuity is fully due to the fluid kinematics. 
Given a fluid distribution with the density p(x\t), one can consider a problem of deter- 
mination of the velocity field u l (x\t) for the fluid motion which preserves the mass of 
its each region. In such a setting the equation of continuity ([H]) is the first order partial 
differential equations to be solved to find the functions of the fluid velocity. This equation 
is not sufficient in general to find 3 unknown functions of the fluid velocity components 
u l {x\€). The dynamics of fluids cannot be therefore reconstructed only on the basis of 
the equation of continuity and the fluid density cannot serve as the physical reason for 
fluid motion. 



9 The Cauchy Equation of Fluid Motion 

The fluid dynamics considers the physical forces, external and internal, to act upon fluid 
particles to cause their motion, and, as a result, a motion of the fluid. The forces acting on 
a moving fluid particle may be of two kinds. Internal, or contact, forces are those acting on 
a fluid particle through its bounding surface, for example, such as fluid pressure. External, 
or body, forces are those acting throughout a medium, for example, such as gravity or 
the electromagnetic force. Due to the Cauchy stress principle |58| over any closed surface 



dH(t) there exists a distribution of a stress vector s l = s 1 ^ ,t;n k ) depending on the 
normal vector n k to the surface, whose action and moment are equivalent to the actual 
internal forces exerted by the fluid outside the region enclosed by this surface upon the 
fluid inside [f^l , [j59f . Such internal forces per unit surface are represented by a 



vector sis = slJx J ,t) directed along an outward normal to the fluid particle surface. 
The external forces per unit mass are represented by a vector F l = F i (x^t) which is 
assumed to be known as functions of position {x 1 } and time t and perhaps also the state 
of the motion of the fluid. In accordance with the Cauchy stress principle the following 
dynamical principle relates the linear momentum P % (t) of an arbitrary fluid region £(t) 



P\t) = j pu l dV (36) 

S(i) 

with the corporate action of the stress vector sV^ and the external forces F l . 

The Principle of conservation of linear momentum The rate of change of linear 
momentum P l {t) ([36| ) of an arbitrary fluid region £(t) is equal to the total resultant force 
acting on the fluid volume within the region 

Jt Pl{t) = Jt J PU%(N= J prdV+ j (37) 

S(t) £(t) dT,(t) 
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By means of the formula ( |53| ) the Principle of conservation of linear momentum can 
be written in the following form: 



P H dV 



P F l dV + 



s\ n) dS. 



(3? 



S(i) 



E(t) 



0£(t) 



In the representation (|3^) an arbitrary moving region E(t) can be replaced, without loss 
of generality, by a fixed region £(to)- An important theorem regarding the stress vector 



s i holds H, B |T 



Theorem 4 (The local equilibrium of the stress forces) For any fluid moving with 
the conservation of linear momentum the stress forces are in the local equilibrium 



lim — 

£(t )-0 D 2 



s\ n) dS = 0. 



(39) 



9£(to) 



Here D is the diameter of the region £(t )- The volume Vs(to) of the region E(t ) is 
Vzito) ~ and the square Sqe^o) of the region surface <9E(t ) Sas^o) — D 2 . 



Proof. This result is derived by division of both sides of ( |38| ) by D 2 and taking limit of 
S(to) - * 0. Then due to the boundedness of both integrands for the volume integrals the 
formula @ gives ( |3|). QED 

Analysis of the structure of the stress vector on the basis of its local equilibrium ( |39"D 
gives the following result p2| , |33l , plfl , |58 |. 



Corollary 3 (The Cauchy formula for the stress vector) The stress vector si-, 
S( n )(a> ? ,t) is a linear function of the normal vector n % = n % {x 3 , t) to the surface <9E(t) 



(40) 



Here the coefficients of the representation form a tensor T u = T^(x k ,t) which is 
called the fluid stress tensor. 

In general the second rank fluid stress tensor T l ^(x k ,t) has no algebraic symmetry in 
its indices [fffil , flST) . 



As a result of the Cauchy formula ([H]) one derives the Cauchy equation of motion 



H, n, p, @ 



Theorem 5 (The Cauchy equation of motion) The Cauchy equation of motion 

du { ■ dT ij 



(41) 



is necessary and sufficient for the linear momentum of any region following the fluid 
particle paths to be conserved. 



16 



Proof. The equation ( [fl| ) follows from the observation that the Cauchy formula for 
the stress tensor (|40| ) and the Gauss theorem enable one to rewrite the Principle of 
conservation of linear momentum (|38D as 



p ^fdv= [ ( P F i + ^-)dv. (12) 



dt J \ dx k 

E(i) E(t) 

Now necessity follows immediately upon substitution of the equation of motion fl4"T| ) into 
(f42|) . Sufficiency can be easily proved from the conservation of linear momentum ( f4~2"[) by 
taking into account that the fluid region E(t) is arbitrary. QED 

The Cauchy equation of motion is the most general equation valid for any fluid motion 
which preserves the linear momentum of its each region regardless of the form of the fluid 
stress tensor. From the point of view of Newtonian mechanics the Cauchy equation of 
motion is the Newton second law of motion stating that the acceleration a l (x\t) ([13] ) of 
a fluid particle with the mass pdV is due to the resulting force of an external force pF l dV 
and an internal force (dT^ /dx^dV acting on the particle. In such a setting the equation 



of motion ([41]) is a system of the first order partial differential equations to be solved to 
find the functions of the components of the fluid velocity u l {x^t). The equation involves 
the density function p(x\t) which satisfies the equation of continuity (|3i~|). Thus, 

The Cauchy equation of motion for the fluid velocity u l 

"{■w +1/ w) = ^ + w (43) 



The equation of continuity for the fluid density p 

| + ^)-0, (44) 



can be considered as a system of four first order partial differential equations for four 
unknown functions of the components of the fluid velocity u % [x^ ,t) and the fluid density 
p(x\t) provided the external forces F l {x\t) and the fluid stress tensor T l ^{x : ,t) are 
given together with the initial and boundary conditions for the velocity and density. The 
dynamics of a fluid is expected to be reconstructed on the basis of solving this system. 
Upon finding the fluid velocity u % {x^ , t), the actual equation motion of the fluid particles, 
either x l = x l (^,t), or = ^(x\t) , is to be found by solving (]T0|). 

The Cauchy equation of motion (f|l]) or (^3] ) can be written in another form by 



introduction of the total acceleration vector A 1 = A l (x 3 } t) which represents the combined 
effect of the inertial force acting on a fluid particle through its acceleration a % = a l (x\t) 
( |T3"D and the effect of the external force F l , 

nil 1 

A* = ^--F l = a l - F\ (45) 
dt 
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The Cauchy equation of motion then reads 



grpij 

PA ' - 7ST < 46 > 

which is transparent from the physical point of view stating that the total acceleration of 
a fluid particle is determined by action of the internal force on the fluid particle. 



10 The Kinematics of a Moving Fluid Particle 



A fluid particle moving in accordance with the equations (|4"3"D and fl4"4] ) is subjected to 
changes in its shape, the value of its volume and it generally rotates with respect to 
its initial, or any other positions {x 1 } at a time t, taken as a reference point. Since it 
is the motion of a fluid particle under influence of the internal and external forces ([43]) 
that affects its state, such changes are due to its changing velocity in dependence of a 
position and a moment of time. Let us consider the tensor dui/dx^ constructed from the 
fluid velocity vector = 5ijV? , which describes the spatial variations in the velocity of a 
fluid particle. This is a second rank tensor has no algebraic symmetry in its indices and 
therefore it can be generally decomposed into the following form: 



duj 
dxi 



(47) 



where the symmetric tensor 9ij = 9ij(x k ,t) called the deformation or rate of strain tensor 
is defined as 



0; 



1 / duj duj 
+ 



2 V dxi dx 

and the antisymmetric tensor Uij = u>ij(x k ,t) called the vorticity tensor is defined as 



(48) 



'J 



1 / dui duj 

2 V dxi dx i 



—Uji. 



(49) 



In many case it is convenient to use the vorticity vector uj 1 

, 1 



(50) 



instead of the vorticity tensor u;^. It is easy to show that uj^uj^ = 0. 

The geometrical meaning of the deformation and the vorticity tensors are given by the 
following theorem [S2"|, ||33|| , due to Cauchy |61| and Stokes |35| . 



Theorem 6 (The Cauchy-Stokes decomposition theorem) An arbitrary instanta- 
neous state of a fluid particle moving along its path may be resolved at each position 
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{^(p)} as a superposition of a uniform velocity of translation u^p) = Ui(x k p ^,t), a di- 
latation d \x l x^9{p)ij\ /dx k along three mutually perpendicular axes and a rigid rotation 



Eijkxiu>f p) of these axes 



Id 1 

Ui(x k , t) = u i{P) + ^g~ [x 3 x k ®(P)jk] + -e ijk x 3 wf F) + 0{x 2 ). (51) 

Here 0(p)ij = %(x^,t) is the deformation tensor at {xl™}, u> k p ~. = oj k p ^(x k p yt) is the 
vorticity vector at {x\ p s\. 

Proof. Let us consider the velocity field u l (x k ,t) of a moving fluid particle in a neigh- 
borhood of its position {x^} at a time t. Since the velocity field is continuous and 
different iable, one has a Taylor expansion of the velocity functions near {X( P \} which 
takes the following form with using the decomposition ([47]) 

Ui(x k , t) = u i{P) + 6 {P)ik x k + uj [P)ik x k + 0(x 2 ). (52) 

The second term in (1527) can be written as 



6{p)ikxk = \h t^w] • ( 53 ) 

This term represents a velocity field normal at each point to the quadric surface x^x k 0^p)jk = 
const which contains the point {x^}. Since a symmetric tensor also possesses three mu- 
tually perpendicular eigenvectors, the eigenvalues of the deformation tensor 9^ measure 
the rates of extension per unit length of a fluid particle at {x| p s} in directions of the 
eigenvectors which can be always taken as three basis vectors for the velocity field (|53|). 



By replacing the vorticity tensor by its vector (|50D the last term in (j52j) it takes the form 

U(P) ik x k = ^e ijk x J uj k P) (54) 

with its immediate interpretation as representing a rigid rotation with an angular velocity 
u\ P) /2. QED 

Thus, the deformation tensor 8ij determines the distortion of the element of fluid 
occupied by a fluid particle along the principal axes defined by the eigenvectors of the 
tensor 6ij, namely a change in the fluid particle's volume and shape. The vorticity tensor 
Uij, or vector a/, determines a rigid rotation of a fluid particle with keeping its volume 
and shape the same. The rotation axis is defined by the direction of the vorticity vector 

There are two important results following from the Cauchy-Stokes theorem. 

Corollary 4 (The rate of change of the fluid path arc) The rate of change in the 
squared element of arc dx l , ds 2 = 5i k dx l dx k , following a fluid particle path is given 

^-(ds 2 ) = 26 ik dx i dx k . (55) 
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Corollary 5 (The rigid motion) For the rigid motion of a fluid when 

9 lJ (x k ,t) = (56) 
everywhere, the fluid velocity u l (x k ,t) is determined everywhere by 

u % = -eijkX^u k + const. (57) 

It is convenient to work sometimes with the trace-free deformation tensor which is 
called the shear tensor cr^ = aij(x k ,t), 

<Tij = Oij - -5ij6, 5 ij aij = 0, (58) 
where 9 = 9(x k ,t) is the scalar of expansion, 

9 = " %■ (59) 

The shear tensor determines the distortion of a fluid particle along the principal axes 
defined by its eigenvalues leaving the volume of the particle unchanged. The magnitude 
of shear a is determined by a =\ a 2 I 1 / 2 , 

^ 2 = ° 2 > °> ( 6 °) 

where a^ 1 = 5 l ^5 kl aik- One can easily prove that a = if and only if = 0. 

The scalar of expansion 9 determines the change of the volume of a fluid particle 
keeping its (spherical) form unchanged. It is useful to define a representative scalar 
R = R(x h , t) of a fluid particle by the equation 

If the expansion scalar vanishes everywhere (|T7|), the fluid is incompressible. If the ex- 
pansion scalar (|59D does not vanish identically for a fluid distribution (|T8|), the fluid is 
compressible. 

The magnitude of vorticity u is determined by to =\ uj 2 I 1 / 2 , 



2 



1 



-u^uJij, uj 2 > 0, (62) 



2 



where u^ 1 = S^S^iJik and it is easy to show that u 2 = uo l uOi. One can prove that u = if 
and only if either cUy = or u l = 0. 
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11 The Conservation of Energy 



On the basis of the Principles of conservation of mass fl28|) and linear momentum (|37|) one 



can establish the transfer equation [32], [51| for the kinetic energy of an arbitrary moving 
fluid region S(t) 

1 



K{t) = - J pSijuVdV. (63) 

E(t) 

Theorem 7 (The energy transfer equation) For any region S(t) following the par- 
ticle paths of the fluid which conserves mass $2$) and linear momentum ( \3l{ ) the rate of 
change of its kinetic energy ( |^5[ ) is equal to the rate at which work is being done on the 
region by external forces reduced by the work which is being done to change the volume 
and the shape of the region and by the work which is being done to rotate the region 



jK>{t)= J p5 ij u i F j dV+ j Sij-u's^dS - J O^T^dV - J uj^T^dV. (64) 

E(t) E(i) 

Proof. The energy transfer equation ( |64"D can be easily derived by using (^), the Cauchy 



equation of motion (|41| ) and a decomposition of the fluid stress tensor T lJ into its symmet- 
ric and antisymmetric parts to produce the last term in accordance with the definitions 
of the fluid deformation tensor 6^- ([47]) and vorticity tensor u>ij ([49|). QED 

The last two terms involving the interaction of the symmetric and antisymmetric parts 
of the fluid stress tensor with the deformation and vorticity tensors can be considered as 
a "dissipation" part. This part of the power in (|64] ) may well be recoverable, but the rest 
must be accounted for as heat. 

Considering now a class of chemically inert fluids without diffusion the total energy 
T(t) of an arbitrary moving fluid region E(t) is defined as a sum of its kinetic energy )C(t) 
( |B"BD and its internal energy £(t), 

T(t) = JC{t)+£(t) = ^ J p5 ij u i u j dV+ J P EdV } (65) 

E(t) E(t) 

where E = E(x l ,t) is the specific internal energy. The thermodynamical interpretation of 
E is different for the compressible and incompressible fluids, (|i~8|) and (|T7|) correspondingly. 
Indeed, for the compressible fluid its density is a function, p(x\t) ^ const, due to the 
continuity equation fl3~T|). Then it is naturally to consider the specific internal energy E 
as a thermodynamical state variable satisfying 

The Postulate of thermodynamical state for the compressible fluids 

TdS = dE + pdr, t = -. (66) 

P 
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Here T is the absolute temperature, S is the specific entropy, p is the pressure, and r is 
the specific fluid volume. Due to this equation of the thermodynamic state (|66| ) the fluid 
pressure p is a thermodynamical variable and there exists in general an equation of state 

p = p(p,S) or p = p(p,T). (67) 

The equation of state is usually chosen or derived on the basis of a phenomenological 
macroscopic and/or a microscopic physical model of the fluid substance. 

For the incompressible fluid p(x p , t) = const due to the continuity equation ( pi]) and 
the pressure cannot be a thermodynamical variable. In this case the hydrodynamic pres- 
sure is not defined as a function of the state of the fluid. Taking one or another definition 
of pressure here, one should remember that it may not be necessarily the physical pres- 
sure measured in experiments. The freedom in choice of the pressure is usually utilized 
to simplify the equations of motion. 

The Postulate of thermodynamical state for the incompressible fluids 

TdS = dE, p = const. (68) 



In general, due to possible internal dissipation processes a fluid can allow a heat 
transfer which is described by a heat flux vector q l = q l (x^ , t) with the physical dimension 
of energy per unit area per unit time. The fundamental postulate of conservation of total 
energy reads as follows |3B|, |5T]] . 

The First law of thermodynamics for fluids (The Postulate of conservation of 
total energy) The rate of change of total energy T(t) ffidj of an arbitrary fluid region 
E(i) is equal to the rate at which work is being done on the region plus the rate at which 
heat is conducted into the region 

j t [K{t) + S{t)} = J P 5 i] u l F 1 dV + j> SijU^dS - j> q i dS l . (69) 

S(t) SE(t) 9E(i) 



The following fundamental theorem establishing the total energy equation now can be 
proved. 

Theorem 8 (The total energy equation) The total energy equation 

pf^+^-g (TO) 

is necessary and sufficient for the First law of thermodynamics for fluids ft6Q) to be satisfied 
for any region which follows the fluid particle paths as that the energy transfer equation 
djfil) for the rate of change of its kinetic energy holds. 
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Proof. Using the definition of the total energy ( |55| ) contained in an arbitrary fluid region 
£(i), Eq. (|35D, the Gauss theorem and the kinetic energy transfer equation (|64[), the First 
law of thermodynamics for fluids (|69| ) takes the form 



E(t) E(t) 

Now the necessity of the total energy equation (|70|) to satisfy Eq. ( |7T|) is evident. The 
sufficiency can be easily seen from (|71]) as a result of arbitrariness of the moving fluid 
region E(t). QED 

The total energy (|70|) in the case of a symmetric fluid stress tensor, = 0, is due 



to Neumann [6^]. Further assumptions on the form of the fluid stress tensor T 13 , the 
heat flux vector q l and the equation of state ( |67D relating the thermodynamical variables 
should be made to apply the First law of thermodynamics for fluids ( |69D for the resolution 
of the energy balance for a particular moving fluid. 

The fundamental postulate of the production of entropy can be formulated for moving 
fluids in the following form g|, B 0- 



The Second law of thermodynamics for fluids (The Postulate of increasing 
entropy) The rate of increase of entropy of an arbitrary fluid region £(£) should not be 
less than the heat conducted into the region divided by the temperature at which the heat 
transfer occurs 

r pSdV <£ L l dS t . (72) 



dt J ~ J T 

E(t) SE(t) 



Theorem 9 (The local increase of entropy) The following inequality for the local in- 
crease of entropy 

P dt = T dx i T 2 dx i 1 ' 



is necessary and sufficient for the Second law of thermodynamics for fluids ( \73j ) to be 
satisfied for any region which follows the fluid particle paths with the conservation of its 
mass 



Proof. Using fl3"o] ) and the Gauss theorem, the Second law of thermodynamics for fluids 
(|72|) takes the form 

dS 1 dq i q { dT\ nr . 

E(t) 

Now the necessity of the inequality (|73|) for the local increase of entropy ([72]) to ensure the 
integral inequality ([74]) follows because an integral of a function everywhere nonnegative 
is nonnegative. The sufficiency can be easily seen from (|74|) as a result of arbitrariness 
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of a moving fluid region £(£), which forces the integrand to be nonnegative everywhere. 
QED 



12 The Fluid Stress Tensor and the Conservation of 
Angular Momentum 

In general, the second rank tensor T l i(x k , t) has no algebraic symmetry in its indices. From 
a physical point of view that means that a fluid can be generally subjected to extraneous 
torques due to action of the external forces F l (x^,t). Such fluids are called polar and 
they are capable of transmitting extraneous stress couples and carrying the body torques 
H, ICO, & If a fluid can only admit the torques arising from the moments of direct 



forces, it is called nonpolar. Thus, certain assumptions concerning the nature of the forces 
exerted across fluid surface elements will specify the properties of the stress tensor. Let 
us accept the following postulate |32|, |6(J about the fluid stress tensor. 



The Boltzmann Postulate The fluid stress tensor T y = T l i(xk,t) is symmetric 



One can prove an important theorem p2| , [[51]], that fluids with a symmetric stress 
tensor T l i(x k ,t) ([F^) conserve the angular momentum L,(t) of an arbitrary fluid region 
E(f) 

Li(t)= J ps ijk x j u k dV. (76) 

E(t) 

The Principle of conservation of angular momentum The rate of change of an- 
gular momentum Li(t) ft7(\ ) of an arbitrary fluid region S(t) is equal to the torque due to 
the total resultant force acting on the fluid volume within the region 



j t Li(f) = j t J p£ijkX j u k dV = J pe ijk x j F k dV+ j 



d 



(77) 



E(t) 



E(t) 



9E(t) 



Theorem 10 (The conservation of angular momentum) For any fluid moving with 
conservation of mass ( \2Q) and linear momentum ( 31_ ) the Boltzmann postulate for the fluid 
stress tensor $7$ ) is the necessary and sufficient condition for the fluid to conserve also 
the angular momentum fiT\ ) for any region following the fluid particle paths. 



24 



Proof. By using fl55| ) and the Cauchy equation of motion flfl] ) one can show that the 
rate of change of the angular momentum ( [76|) can be expressed in the following form: 

^Li(t) = y pe ijk x J u k dV = J pe ijk x J ^_u k dV = 

E(t) E(t) 

y pe ijk x j F k dV + y - y e ijk T jk dV. (78) 

s(t) as(t) s(t) 

Since Eij k T^ k {x l ,t) = if the fluid stress tensor is symmetric, the Boltzmann postulate ( [75]) 
necessarily brings the Principle of conservation of angular momentum (|77|). Sufficiency of 
the Boltzmann postulate follows from (|77f) and (|78|) due to arbitrariness of the fluid region 
E(i) when the vanishing of the last integral in ( |78|) requires £ij}.T^ k {x^t) = everywhere. 
QED 

The interpretation of the components of a symmetric fluid stress tensor is well-known 
(see for example, [RH], pT|). 



13 The Navier-Stokes Equation and Its Physical Sta- 
tus 

The different physical hypotheses about the nature of moving fluids and forces causing 
their motion bring different analytical forms for the fluid stress tensor. As a result, 
historically various forms of the Cauchy equation of motion ( f43|) have been discovered to 
describe the dynamics of various fluids. By restriction to the class of fluids without heat 
transfer, 

q i {x j ,t) = 0, T = const, (79) 
there are three forms of the equations of motion of interest and relevance here. 
The Euler equation of motion for the perfect fluid T u = — S^p 



p2| , [51 1, [p4 ] . Here p = p(x J ,t) is the fluid pressure which depends on the thermody- 



namical state of fluid. It can be positive, p > 0, which tends to compress a fluid particle, 
and negative, p < 0, which results in an expansion of a fluid particle. To close the sys- 
tem of the four first order partial differential equations (pH|) and (^Up for four unknown 
functions of the components of the fluid velocity u l [x\ t) and the density p(x l , t) provided 
the external forces F % (x^,t) and the initial and boundary conditions for the velocity and 
density are given, one must specify in addition to the system the of state ( |67D which does 
not depend on temperature due to (|79|) 



25 



The equation of state 

p = p(p) or p = p(p). 



(81) 



The Stokes equation of motion for the Stokean fluid T' lj = (—p + a)5' lj + j39 l ° + 



p2| , p5| , |5l| . Here a = a(x\t), (3 = f3{x\t) and 7 = r y(x i ,t) are the scalar functions 
which depend on the thermodynamical state of fluid as well as on the principal invariants 
of the symmetric deformation tensor l] = 9 lJ (x k ,t), 6 lj = 5 lk 5 jl 6 k i ([4|). They describe 
definite types of viscous response of the fluid. For example, a choice of p, a, [3 and 7 so that 
T*- 7 is linear in 6^ provides the classical Cauchy-Poisson law of viscosity. The Stokean 
fluid is the most general type of a nonelastic fluid. To close the system of Eqs. (|44]), 
(|8~TD and fl82|) for four unknown functions of the components of the fluid velocity ■u i (x jf , t) 
and density p{x\t) provided the external forces F l (xi,t) and the initial and boundary 
conditions for the velocity and the density are given, one must specify in addition the 
viscosity functions a(x\t), (3(x l ,t) and r y(x\t). 

The Navier-Stokes equation of motion of classical hydrodynamics for the New- 
tonian fluid T ij = 5 ij (X9 -p) + 2[M6 ij 



P3j], P5] , [ fxTf . Here A = X(x\t) and [i = [i(x\t) are the scalar functions of the 
first and second coefficients of viscosity. For the case of constant viscosity coefficient, 
X(x\t) = const and p,(x\t) = const, the Navier-Stokes equations becomes much simpler 

\ ot ox 3 J ox 3 oxWx K ax K ax l 



Its further simplification for the case of incompressible fluid (|T7|) has been especially 
thoroughly studied in connection with the phenomenon of fluid and gas turbulence [H| , 



The Newtonian fluid is a linear Stokean fluid with the Cauchy-Poisson law of viscosity. 
For a Newtonian fluid its stress tensor depends linearly on the rate of fluid deformation. 
It is generally believed that the Newtonian fluid is an adequate physical model when the 
fluid deformations are relatively small in comparison with, say, the ratio of some reference 
speed and reference length. However, this physical hypothesis is not well understood as yet 
. It is not to be derived from experiments, nor can it be proved by abstract modelling. 



The Navier-Stokes equation (|83| ) is the dynamical equation of classical hydrodynamics. 
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Since their discovery they successfully describe various regimes of fluid motion including 
turbulence for a wide class of physical gaseous and liquid substances. 

A closed system of the equations of classical hydrodynamics without heat transfer 
includes Eqs. ([44]), ( |81~D and (|84]) for four unknown functions of the components of the fluid 
velocity u l (x^,t) and the fluid density p(x\t) provided the external forces F l (x^,t), the 
viscosity coefficient functions A(x l ,t) and fj,(x l ,t) and the initial and boundary conditions 
for the velocity and the density are given. 

The mathematical difficulties in solving the Navier-Stokes equation are well known, 
see [^TJ for a discussion and further references. Because the 3-dimensional Navier-Stokes 
equation is not a well-posed partial differential system, no theorem of existence of smooth, 
physically reasonable solutions is not known as yet. These problems originate in the 
nonlinear structure of the Euler equation fl80f ) and the Navier-Stokes equation fl84]) due 
to the presence of the so-called inertial term u^du 1 /dx J ' . These equations are among the 
very few equations of mathematical physics for which the nonlinearity arises not from 
the physical attributes of the system, such as a physical interaction, but rather from the 
kinematical aspect of fluid motion. Another source of the difficulties for the Navier-Stokes 
equation is due to the viscosity term which brings the second order operator into the first 
order equations, which involves the singular perturbations. 

As far as the physical status of the Navier-Stokes equation is concerned it must be 
regarded as the macroscopic equation to be able to describe the dynamics of fluids con- 
sidered as continuous physical substances. The flow of fluids at the microscopic level is 
governed by the statistical mechanics of fluids. The appropriate statistics is given by 
the solution of the Boltzmann equation which represents the evolution of the distribution 
function in dependence on positions and velocities of colliding particles as a result of 
thermal excitation. The collisions are described by the so-called collision integral which, 
in principle, contains all information about interaction between particles and it forms the 
collision operator. One must evaluate this operator to obtain a set of macroscopic equa- 
tions for one or another physical configuration of interacting (colliding) particles. Only 
for the case of dilute gases when one is limited to the case of the evolution of a single 
particle distribution with only binary collisions taken into account it has been possible to 
find a satisfactory solution given by Chapman |k| and Enskog fl6l| . It has been shown 
that the Boltzmann equation with this collision operator gives the equation of the con- 
servation of mass and the equation of the conservation of linear momentum which is the 
Navier-Stokes equation for an incompressible fluid []65[| . Although such a derivation has 
been carried for dilute gases, a derivation for liquid fluids when the binary collisions are 
almost of negligible significance and all liquid physics is due to collisions between particle 
clusters is still an open problem. In this situation the Navier-Stokes equation (34) is taken 
as a reasonable model for liquid flows at least for the simple constant viscosity. 
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